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^^ • Abstract. We consider an elastic system containing a small region where the 

^vi ' density is very much higher then elsewhere. Such system possesses two types 
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expansions of global eigenvibrations for ordinary differential operator of the 
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Introduction. Heterogeneous systems give rise to new effects that do not reside 
in the separated system parts. For instance, problems with local density perturbation 
are characterized by a presence of local and global vibrations. For the first time the 
effects have been described by E. Sanchez-Palencia [6,7]. Vibrating systems with a local 
V") , mass perturbation are investigated starting from 0. A. Oleinik [3,4]. The effect of local 

vibrations has been studied for wide class of the systems, for instance see [1,2]. The global 
vibrations have been remained to be weakly analyzed so far. A complete asymptotic 
expansions of the global eigenvibrations for one-dimensional system of the forth order 
with locally perturbed density are constructed in this work. 

O ' 

qq ' 1. Problem statement. Let a differential expression L be given by 

^^ ' Lv = (ko{x)v")" — (ki(x)v')' + k2(x)v, 

where the functions fco > 0, fci, &2 > are smooth at [a,b]. Denote by [ifja—c a jump of 
k>( " a function v at point c. We investigate asymptotic behaviour as e — » of eigenvalues A e 

^N , and eigenfunctions u e of the problem 

_C3 J Lu e - X E p(x)u e = 0, x € (a, -e)U (e, 6), (1) 

Lu E - \ E e~ m q(x / e)u e = 0, x€(— e,e), (2) 

u e (a) = Us (a) = 0, u e (b) = u' e (b) = 0, (3) 

[Ue] x=±e = [t4]«=± e = [<%=± S = K']«=±« = °- ( 4 ) 

For each fixed e > the problem possesses a countable set of eigenvalues. Bihaviour 
as e — > for each eigenvalue A e and corresponding eigenfunctions u E depending on a 
value of real parameter m is investigated in [1]. In the case m > 4 problem (1) - (4) 
possesses local eigenvibrations 1 with corresponding eigenvalues A^, = 0(e m_4 ), e — » 0, 
and with eigenfunctions u E n , which localize itself in the interval of density perturbation 
(— e, e), rapidly vanishing outside the interval. Nevertheless, qualitative behaviour of the 
vibrating system is not yet described completely by local eigenvibrations. Models with 
concentrated masses possess also global vibrations [7]. As indicated below, the global 



x We use a term "eigenvibration" to denote a pare of eigenvalue and corresponding eigenfunction 
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vibrations are supported by eigenfunction sequences with nontrivial limits u e n r e \ — > vo for 
eigenvalues A^( s -) — » Aq > with n(e) — > oo. Dependence n(e) is a discrete one that causes 
construction of asymptotics along certain sequences of the small parameter only. A family 
of the sequences is bound by a deformation parameter, which is present in the lower terms 
of the constructed expansions. Asymptotics depend on the value of m. We choose m — 8 
as a pattern case. 

2. Asymptotics of global vibrations: the leading terms. We seek the asymptotic 
expansions of the eigenvalues A e and the eigenfunctions u e of problem (1) - (4) in the form 

oo 

A e ~^>%, (5) 

i=0 

oo 

Ue{x) ~ 2je ! Wi(s), a; G (a, — e) U (e, b), vo ^ 0. (6) 

i=0 

In order to explore the function u e in the region (— e,e), we consider problem (1) - (4) in 
variables £ = e~ 1 x. Taking into account 



fc i (eO = £(^) J '*i , ' ) (0)(i! 



j=o 



we transform the differential expression I to I E = e 4 \^ e J C(j), where 

c( ~ = A^(0)rf^_tf _ k { r 2) (0) d j_ 2 d fc^~ 4) (0) ,--4 
UJ j! dC 2? rfC 2 (i-2)! dr ^ (J -4)! ? ■ 

We use notation k„ (x) for the i-th derivative of a function k n (x) in case if i £ N, k„ (x) — 
k n (x) and k { n j) (x) = if j < 0. 

Hence, the eigenfunction u e , which corresponds to A e , is a solution of the equation with 
small parameter nearby the highest derivative 

c 8 L e {7 £ - X e q(OUs =0, U e = 17.(0, £ G (-1, !)■ (7) 

According to the method of WKB-approximations [4], we seek the solution of (7) as linear 



combination of the series e £ "' \e 8 <Zi(£)- Equality (7) is guarantied by the choice of 

i=0 

functions S and at- In particular, the phase function S is a solution of the eikonal equation 
k (0)S' 4 - X q(0 = 0, ££(-1,1). (8) 

Therefore we fix 

S(C) = (Aofc " 1 (0)) 1/4 |? 1/4 W df 
-l 
and, without loss of generality, we derive the eigenfunction u £ in the form 

oo 

u e (eO= f U.(£) = e*J2s i (f i {t),N{t, e - 1 S)), {€(-1,1), (9) 

i=0 

with fi being a vector-function with values in K 4 , (•, •) being the standard scalar product 
in R 4 , and the operator N mapping [— 1, 1] x C°°[—l, 1] in R 4 according to 

JVtt,T) = (oo8T(0,sinT(0,exp(-T(0+T(-l)),exp(T(0-T(l))). 

The choice of the multiplier e 4 and the argument shifts in exponents of series (9) is related 
to a certain normalization of the eigenfunction u e , which is described in Section 4 (all other 
eigenfunctions differ by a constant multiplier). 
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We construct formal asymptotic expansions (5), (6), (9) satisfying all conditions of 
problem (f) - (4). In particular, equality (I) holds if 

Lv - \op{x)v = 0, x G (a, 0) U (0, b), (10) 

i 

Lvi — \op(x)vj = p(x) y^ \jVi-j, x 6 (a, 0) U (0, 6). (11) 

j=i 
Let consider the action of the differential expression L e on the function U e . We note that 
— N(£,t) = t'(£)TN(£,t), where the matrix 

at; 

T = ( T 2 )> Tl = { 1 ) ' T2= { 1 
has the property T* = T 3 , and T 4 is a unit matrix. Counting 

we get 

oo i 4 

L e Ue = E £l " 8 E(E^(j -o/i-i(0,^(e,e -1 s)), 

»=0 3=0 (=0 

where £„(fc) are the differential expressions of the n-th order with coefficients depending 
on S" and £ fc . In particular, £ (i) = /^(OX* 1 )"^' 4 ^, and 

We recall that the differential expressions C„(t), n = 0, . . . , 4, are equal to zero one for 
t < 0. Hence, equation (2) yields 

(£o(0) - Xoq)fi + (£o(l) +£i(0) - Aig)/<-i = Xi, 

i 4 

with Xi = — /^ /(AO — ~~ ^jl)fi-j f° r * > 2, and xt = for £ < 2. The expression 

3=2 (=0 

£o(0) — Ao<? is equal to zero because it is the left-hand side of the eikonal equation given 
by (8). Therefore the coefficients /, of expansion (9) solve the system of the first-order 
differential equations 

OCo(l)+£i(0)-Aig)/t = X *+i, i = 0,1,2,.... (12) 

Boundary conditions (3) give 

Vi{a) = v'i(a) = 0, Vi(b) = v'i(b) = 0, i = 0,1,2,.... (13) 

Interfacial conditions (4) applied to (6) and (9) yield 



\ at at at at / 



E^WrM^io) = </ ! - 4 (±l),iV(±l, e - 1 S)>, 



3=0 

E(±l) J (j!)" 1 ^ 1) ( ± °) = <(S'T 3 / ! -2 + /;- 3 )(±l),iV(±M~ 1 S)), 

3 = 

i 

E(±1) 3 0'!) _1 ^ 2) (±0) = {(S' 2 T 2 f i +D i )(±l),N(±l,e- 1 S)), 

Fi(±0) = ((^'T/i + E i )(±l),N(±l,e- 1 S)), 
where we used notation 

fj=0 for j<0; 

A = 2S'T 3 fU + 5"T 3 /i-i + .A'-2 for i = 0, 1, . . . ; 



(14) 



(15) 
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Fi take zero values for I = 0,1, and 

F i (±0) = X)(±l) i (j!)" 1 «H.'Ji ) a (±0) for i = 2,3,...; 

and ^ - { S ' 2 T, + S T, S ' + T, S ' 2 ) T ^ + (*'$ + |*'| + P' 

T' A f 3 -2 + f"U for i = 0, 1, . . . . 

According to conditions (10), (13), (14) for j = 0, the leading terms Ao and Vq of 
expansions (5), (6) are correspondingly an eigenvalue and eigenfunction of the problem 

Lv - X p(x)v = 0, x G (a, 0) U (0, b), . 

vo(a)=v' (a) = 0, v (0) = v' (0) = 0, v (6) = « (&) = 0. l J 

We restrict ourself to considering only simple eigenvalues of three-point Dirichlet prob- 
lem (16). Nevertheless, there are situations when all the eigenvalues has multiplicity more 
then 1 (for instance, a — —b and p is an even function). We fix a simple eigenvalue Ao 
of (16) and the corresponding eigenfunction vq such that Vq(x) = for x € (0, b) and 

PVq dx — 1. 

Equalities (11), (13), (14) for i = 1 give 

Lvi — \op(x)vi = Xip(x)vo, x £ (a,0) U (0,6), 
«i(o)=vi(o) = 0, wi(-0) = 0, «i(-0) = «J'(-0), (17) 

v 1 (+0) = v' a (+0)=0, vi(b)=v[(b) = 0. 

Three-point problem (17) has solution if and only if the parameter Ai takes value 

Ai = fco(0)«o'(-0) 2 . 

We fix a solution V\ of (17) such that / pv\Vo dx = 0. One can note that vi(x) = for 

xe (o,fo). 

Conditions (12), (15) for i = yield that the leading term /o of series (9) satisfies the 
problem 

fo = MZ)fo, £e(-i,i), (18) 

(T fl /o(±l),JV(±l,e- 1 S)> = v*, (Tfo{±l),N(±l,e- 1 S)) = 0, (19) 

with <r + = 0, a~ = S'(— 1)~ 2 v'q(— 0). The matrix of the linear homogeneous system of 
the first-order differential equations (18) is 

with 

77 = -3g'(C)(8g(C))- 1 , 9 = 1/4 (A ~ 3 fc " 5 (0)g(C)) 1/4 (Aifco(O) - \ o k' Q (0)O- 

On the one hand, problem (18), (19) depends on the small parameter, on another hand, 
it is a boundary value problem, and hence, it is an ill-posed one. Both of the difficulties can 
be solved by exploring problem (18), (19) along with similar problems, which appear below, 
along discrete sequences of the small parameter e; — ► 0. Then "along" these sequences 
the problems have a unique solution, which is independent of ei up to the exponentially 
small terms. In particular, this property is consistant with a discrete character of the 
high-frequency-vibrations effect, which is under consideration. 
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Proposition 1. Let w : [— 1, 1] — > R 4 be a smooth vector-function, and a be a vector in 
M 4 . There exists a small parameter sequence {ei}fZi such that the problem 

y'(Z,e) = A(Oy(t,e)+w(0, Ce(-M), (20) 

(T 2 y(-l,e), N(-l,e- 1 S))=a 1 , (Ty(-l,e), Ni-l^S)) = a 2 , 

(T 2 y(l,s),N(l,s- 1 S))=a 3 , (Ty(l,e), N(l, e^S)) = a 4 , ^> 

has a unique solution y(-, ei) for each I = 1, 2, . . . . 

The family of solutions {y(-,ei)}fZi satisfies the estimate 

M-^^-y^Wc^Ce^p 1 " 1 

with y* being a smooth vector function on [— 1, 1], and with positive constants C and M 
that are independent of e. 

Proof. System (20) has a fundamental matrix 

* = „-*/*( ^l \ 9l= ( COSQ SillQ 

V $2 / ' V — sin a cos a 

- ( e aC_1)_a 

2 ~ V e a ~ a 

with a function 

«(0 = i(A - 3 fco" 5 (0)) 1/4 |9 1/4 ( i )( Alfc o(°) - Ao*o(0)t)«. 

-i 

Then the general solution of (20) can be written as 

i/(O = *(0G8 + M0). 

where j3 is a constant vector, and 



MO = / $ -1 (*M*) d *- 



Let further 

y(£,e) = *«)(& + MO) 

be a solution of boundary value problem (20), (21). It is easy to check that 

**(0^(€,e _1 5) =<T 3/8 (£)^(£,7 e ), where 7e (£) = e" 1 ^) + a(£), 
and $ f is transposed to $. Therefore we get 

(y^,e),N^,e^S)) ^ q' 3/8 (0{l3. + h(0,N^, lE )), 

that allows to rewrite boundary conditions (21) to the form 

{j3e,T 2 N{-l, le )) =m-cn, 
{p e ,T 3 N{-l, le )) =m-a 2 , 



(/3 £ ,T 2 7V(1, 7£ )) =m+a 3 -{h(l),T 2 N(l,^))), 
</3 £ ,T 3 iV(l, 7£ )) = m + «T4 - (h(l), r 3 iV(l, 7£ ))), 

where m ± = g 3//8 (±l). Here we have used the equalities T$* = $'T and ft(— 1) = 
Therefore, the vector /3 e has to be a solution to a linear algebraic system with matrix 



(22) 



G(7*(l)) = 



/ -1 1 e">M \ 

' -1-1 e-^ 1 * I 



-C0S7 e (l) -sin 7e (l) e -' s{1 > 1 

V Sin 7s (l) -cos 7e (l) -e"T £(1) 1 / 

Note that the determinant 

detG( 7e (l)) = -2cos 7e (l) + 2e" 7e(1) (2 - e" 7e(1) cos 7e (l)) 



-1 





1 


°\ 





-1 


-1 





- cos 6 


— sin 8 





1 


sin 8 


— cos 8 





1 / 
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becomes zero on an infinitely small sequence e (but not for any e > since 7 e (l) -^ooas 

£ -►(}). 

We fix a number 8 from the interval [0, 2tt) such that 8 is different from 7r/2 and 37r/2. 
We construct the sequence ei using the set of conditions 7e;(l) = 8 + 2nl for I — 1, 2, . . . . 
In other words, 

e 'ffl = rf + 2%-a(l) < 23 > 

for all I > lo, where lo is the smallest natural number such that the denominator of 
fraction (23) is positive. We use notation 7; = r ye l . Using 7; — » 00 as ei — > 0, we can treat 
the matrix G(yi) as an exponentially small perturbation of matrix 

/ 

G& = 

V 

In the same manner iV(l,7;) = ( cos5, sin<5, e - ^' 1 ', 1 ), then the right hand side of sys- 
tem (22) also is an exponentially small perturbation of the vector 

g = (m~ai, m~~<T2 : m + a s - (h(l), T 2 Ni), m + a 4 - (h(l), T 3 Ni}), 

where Ni differs from N(l,~/i) only by the third coordinate being zero. 

For the chosen value S the matrix Gs is a non-degenerate one. Then according to the 
results of perturbation theory in a finite dimensional space, we have 

||/3 £i -/3*|| K 4<C7e- wl \ 

where /3* is a solution to the system GsP = g. Let J/*(£) = C &(C)(/^* + M£))i then 

ll»(-,eO-tf*lloi<ll*lloi||^,-ja.|lB* ) 

with ll^llc 1 = max (II^COII + ll^'(OII)- To finish the proof, it remains to note that 

£€[-M] 
7e(1) > Me^ 1 with a certain positive constant M. D 

Remark. The function y* is referred to as the principal solution of problem (20), (21) 
since constructing asymptotic expansions in power scale of e we can neglect exponentially 
small terms. Nevertheless, the choice of the sequence ei, and then the choice of the 
solution y„, is not unique since it depends on 8. As we mentioned at the beginning, 
the latter is connected to the fact that the discrete approximation of global vibrations 
allows deformation, hence the approximation is not unique. A presence of the deformation 
parameter 8 in asymptotics corresponds to the problem content. 

We come back to exploring problem (18), (19), which is problem (20), (21) with the 
right-hand side w — and a — (S'(— l)~ 2 v'o(— 0),0, 0, 0). The system is homogeneous, 
and then /o(£) = < I > (5)/3o, where 

Po = 1/2 g 3/8 (-l)S'(-l)- 2 ^(-0) (tg<5 - 1,-tgS- l,tg<5+ 1, -1/ cos 8) 

is a solution of the corresponding linear system with matrix Gs- 

Hence, we have found the terms Ao, vo, Ai, Vi, /o of expansions (5), (6), (9). We 
recall that /o depends on value 8. The construction of complete asymptotic expansions 
for solutions to problem (1) - (4) is conducted on the sequence £$ = {si(5)}f° that is 
chosen accordingly to (23). 

3. Complete asymptotics of global vibrations. Let find coefficients Ai, Vi, fi-i 
of asymptotic expansions (5), (6), (9) for i > 2. Using conditions (11), (13), (14) we 
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construct a boundary value problem for the general term Vi in the form 

i 

Lvi — \ Q p(x)vi — p(x) 22 ^3 V i-3> x £ ( a i 0) u (0, b), 

j=i (24) 

v t (a) = v'i(a) = 0, vi (6) = v[ (b) = 0, 
«i(±0) = Vi(±0), »<(±0) = Wi(±0). 

The right-hand sides of the boundary conditions in problem (24) with the precision up to 
exponentially small terms are 

i-l 

V t (±0) = («- 3/8 $- 1 /i-4(±l),iV ±1 )-^(±l) i (i!)- 1 V |i ) j (±0), 

j=o 

and 

i-l 

Wi(±0) = {q-^Q-'iS'T 3 /^ + fi_ 3 )(±l),N ±1 ) - ^(±l) J O , 0" 1 «H-i 1) (±O), 

3=0 

where JV_i = (1, 0, 1, 0), and the vector iVi was introduces during the proof of Proposition 
1. Since Ao is an eigenvalue of problem (16), boundary value problem (24) could have no 
solution. The necessary and sufficient condition for existence of the solution is equality 
Xi = (kov'dWi — (koVo)'Vi)(—0). We fix the solution vt to problem (24) such that 

o 

v^o — 0. 



Then we find fi-i from the problem 

fl = A(t)ft + 1/4 fc5" 1 (o)5'- a r 8 xifi ) £ e (-1, i), 

(r 2 /«(±l), N(±l,e- 1 S)) = -((g- 3 / 8 S'- 2 $- 1 A)(±l),iV ±1 )+ 

i 

S'(±l)- 2 ^(±ir(i!)- 1 ^(±0), (25) 

<T/ l (±l),iV(±l,e- 1 S)) = 

S'(±l)- 3 (f,(±0) - ((g-'/'S-^JCilJ.JVii)) . 

According to Proposition 1, the right-hand side of the system is a smooth function. Then 
for the sequence of small parameter £g there is a solution to problem (25). 

Hence the algorithm of constructing coefficients of expansions (5), (6), (9) and a se- 
quence of small parameter is 

Ao — » wo — ► Ai — > vi — > £ s — » /o —►•••—► Xi —* Vi —> fi-i .... 

We recall that all the coefficients, except Ao, Vo and Ai, Vi, depend on the value of 
parameter S. In other words, to each S £ [0, 27r), S ^ n/2, 3ir/2 we assign the corresponding 
asymptotic series that approximate the same form of global vibrations vo . 

4. Justification of asymptotics. The idea of justification is, using asymptotic 
series (5), to choose exactly the eigenfunction sequence of problem (1) - (4) such that the 
sequence models global vibrations and is approximated by series (6), (9). 

We fix numbers n G N and 8 6 [0, 2-n) different from 7r/2 and 3-7r/2. We introduce 
a number sequence {AJ™ } £ c p, and a sequence of functions {ui } £ c p.. Namely, for 
each element e of the discrete set £$ we assume 

A] n} = Ao + eAi + . . . e n X n , 
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!vo(x) +svi(x) H +e n v n (x), x e (a, — e) U (e, 6), 

n+2 
e 4 ^eV/iCe-HJVCe- 1 !^- 1 ^), z G (-e,e), 
i=0 

where the values Ai, functions Vi, vectors ft and the set £s are constructed in Sections 2 
and 3. 

Proposition 2. There exists an eigenvalue sequence {A e } £ ^ p such that 

\X { s n} -X E \<C n e n+1 , ee£s. 

Proof. Denote by p E the density of the original problem, which equals e~ s q{x/e) 
on the interval (— e,e), and is p(x) outside the interval. Problem (1) - (4) considered in 
the weight space Lz^pe, (a, &)) is equivalent to the problem for the self-adjoint operator 
C e = P7 X L: 

C e u E — X E u e = 0. 

The domain of the operator C E is 

D(C e ) = {</?£ H\a, b) ■ <p(a) = <p'(a) = 0, <p(b) = <//(b) = 0}. 

The spectrum of the operator £, e is discrete [1]. 

The function i4 is not an element of the space D(C e ) since it has discontinuities 
at the points x — ±e. Nevertheless, we state existence of the function </4 such that 
u\ n + tpe G D(Ce), moreover </4 is equal to zero on the interval (— e,e), and 



max N^ 

x£(a,b) *— ' 



4 <f 



dx i 



< C n e n+1 . (26) 



i=0 

We introduce function 

V e {n} = M*4 n} + v4 n} ), 

where the factor k £ is chosen accordingly to 

\\ V e " \\L 2 (p e: (a,b)) - 1. 

We remark that the value n e is separated from zero by a positive constant that does not 
depend on e. 

Taking into account (26), we can show that 

\\r \A n } \ {"■}T/-*f ri } II <r k C r n + 1 - e p r 

II IIL 2 (p e ,(a,6)) 

where the constant T^n is independent of e. Therefore there exists an eigenvalue A e of the 
operator £ £ such that 

\r -\ { E n} \<K n E n+1 . u 

Using asymptotics of each eigenvalue A^ we can prove that in a certain de 4 -vicinity of 
points A e , which were chosen accordingly to Proposal 2, we don't have other eigenvalues of 
the initial problem. As consequence, for n > 5 we have only one eigenvalue A e that satisfies 
Proposition 2. We choose the sequence of exactly these {A e } £ ^ p and the corresponding 
normalized eigenfunctions {u e } £ p c . Then 

\\ue-KeV { e n} \\ L2{pe){ajb)) < K n d-\ n -\ ll««lli 3 (p e ,(a ) 6))= 1 - 
This inequality complete the proof of 
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Theorem. On the sequence £s the chosen functions u E satisfy the estimates 



Me (a;) - k £ y^^e k v k (x] 



<C n e 



n+l 



L 2 (fi e ) 



where region Q s states for (a, —s) or (s,b), and the estimate 



u £ {e£) - £ 4 « s J>* </*(£), N(£, e^S)) 



<C n e 



n+l 



i 2 (-l,l) 



for n = 0, 1, . . . . The coefficients Vk, fk are the same as at the beginning of Section 4- 
The normalizing multiplier k £ tends to 1 as e — » 0. 
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